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^ ; Abstract 

H ' The spin superfluidity - superfluidity in the magnetic subsystem of a condensed matter - is man- 

ifested as the spontaneous phase-coherent precession of spins first discovered in 1984 in He-B. This 
superfluid current of spins - spin supercurrent - is one more representative of superfluid currents 
\ known or discussed in other systems, such as the superfluid current of mass and atoms in superfluid 

^He; superfluid current of electric charge in superconductors; superfluid current of hypercharge in 
' 5^ '■ Standard Model of particle physics; superfluid baryonic current and current of chiral charge in 

^ ■ quark matter; etc. Spin superfluidity can be described in terms of the Bose condensation of spin 

waves - magnons. We discuss different phases of magnon superfluidity, including those in magnetic 
trap; and signatures of magnons superfluidity: (i) spin supercurrent, which transports the magne- 
I tization on a macroscopic distance more than 1 cm long; (ii) spin current Josephson effect which 

shows interference between two condensates; (iii) spin current vortex ~ a topological defect which 
■ is an analog of a quantized vortex in superfluids, of an Abrikosov vortex in superconductors, and 

^ . cosmic strings in relativistic theories; (iv) Goldstone modes related to the broken U{1) symmetry 

Q ! - phonons in the spin-superfluid magnon gas; etc. We also touch the topic of spin supercurrent in 

general including spin Hall and intrinsic quantum spin Hall effects. 

K*" \ PACS numbers: 

' Keywords: Bose-Einstein condensation, magnon, spin supercurrent, superfluid '^He, Q-ball, spin-Hall effect 
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I. INTRODUCTION 



Nature knows different types of ordered states. 

One major class is represented by equilibrium macroscopic ordered states exhibiting spon- 
taneous breaking of symmetry. This class contains crystals; nematic, cholesteric and other 
liquid crystals; different types of ordered magnets (antiferromagnets, ferromagnets, etc.); 
superfluids, superconductors and Bose condensates; all types of Higgs fields in high energy 
physics; etc. The important subclasses of this class contain systems with macroscopic quan- 
tum coherence exhibiting off-diagonal long-range order (ODLRO), and/or nondissipative 
superfluid currents (mass current, spin current, electric current, hypercharge current, etc.). 
The class of ordered systems is characterized by rigidity, stable gradients of order parame- 
ter (non-dissipative currents in quantum coherent systems), and topologically stable defects 
(vortices, sohtons, cosmic strings, monopoles, etc.). 

A second large class is presented by dynamical systems out of equilibrium. Ordered states 
may emerge under external flux of energy. Examples are the coherent emission from lasers; 
water flow in a draining bathtub; pattern formation in dissipative systems; etc. 

Some of the latter dynamic systems can be close to stationary equilibrium systems of 
the flrst class. For example, ultra-cold gases in optical traps are not fully equilibrium states 
since the number of atoms in the trap is not conserved, and thus the steady state requires 
pumping. However, if the decay is small then the system is close to an equilibrium Bose 
condensate, and experiences all the corresponding superfluid properties. 



A. BEC of quasiparticles 

Bose-Einstein condensation (BEC) of quasiparticles whose number is not conserved is 
presently one of the debated phenomena of condensed matter physics. In thermal equilibrium 
the chemical potential of excitations vanishes and, as a result, their condensate does not 
form. The only way to overcome this situation is to create a dynamic steady state with a 
conserved number of excitations as a non-equilibrium system. Thus the Bose condensation 
of quasiparticles belongs to the phenomenon of second class, when the emerging steady state 
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of the system is not in a full thermodynamic equilibrium. It is decaying, but the loss of 
quasiparticles owing to their decay can be compensated by pumping of energy. 

Formally BEC requires conservation of charge or particle number. However, condensa- 
tion can still be extended to systems with weakly violated conservation. For sufficiently 
long-lived quasiparticles their distribution may be close to the thermodynamic equilibrium 
with a well defined finite chemical potential, which follows from the quasi-conservation of 
number of quasiparticles, and the Bose condensation becomes possible. Several examples of 
Bose condensation of quasiparticles have been observed or suggested, including phonons 
excitons [2], exciton-polaritons js], H]- The BEC of quasi-equilibrium magnons - spin waves 
- in ferromagnets has been discussed in Ref. ^ and investigated in jl, [^. 

The coherent spin precession observed in superfiuid phases ^He also can be described in 
terms of magnon BEC. The coherent spin precession was discovered first in ^He-B, where it 
is was called the Homogeneously Precessing Domain (HPD) js], llf. This is the spontaneously 
emerging steady state of precession, which preserves the phase coherence across the whole 
sample even in the absence of energy pumping. The crucial property of HPD is that pre- 
cession acquires a coherent phase throughout the whole sample even in an inhomogeneous 
external magnetic field. This is equivalent to the appearance of a coherent superfiuid Bose 
condensate. 

In the absence of energy pumping this HPD state slowly decays, but during the decay 
the system remains in the coherent state of BEC: the volume of the Bose condensate (the 
volume of HPD) gradually decreases with time without violation of the observed properties 
of the spin-superfiuid phase-coherent state. A steady state of phase- coherent precession can 
be supported by pumping. But the pumping need not be coherent - it can be chaotic: 
the system chooses its own (eigen) frequency of coherent precession, which emphasizes the 
spontaneous emergence of coherence from chaos. 

HPD is very close to the thermodynamic equilibrium of the magnon Bose condensate 
and exhibits all the superfiuid properties which follow from the off-diagonal long-range order 
(ODLRO) of the coherent precession. After discovery of HPD, several other states of magnon 
BEC have been observed in superfiuid phases of ^He, which we discuss in this review, 
including finite magnon BEC states in magnetic traps. 



B. Spin superfluidity vs superfluidity of mass and charge 

Last decade was marked by the fundamental studies of mesoscopic quantum states of 
dilute ultra cold atomic gases in the regime where the de Broglie wavelength of the atoms is 
comparable with their spacing, giving rise to the phenomenon of Bose-Einstein condensation 
(see reviews (lol.[ll|). The formation of the Bose-Einstein condensate (BEC) - accumulation 
of the macroscopic number of particles in the lowest energy state - was predicted by Einstein 
m 1925 [l2|. In ideal gas, all atoms are in the lowest energy state in the zero temperature 



limit. In dilute atomic gases, weak interactions between atoms produces a small fraction of 
the non-condensed atoms. 

In the only known bosonic liquid ^He which remains liquid at zero temperature, the BEC 
is strongly modified by interactions. The depletion of the condensate due to interactions 
is very strong: in the limit of zero temperature only about 10% of particles occupy the 
state with zero momentum. Nevertheless, BEC still remains the key mechanism for the 
phenomenon of superfiuidity in liquid ^He: due to BEC the whole liquid (100% of "^He 
atoms) forms a coherent quantum state at T = and participates in the non-dissipative 
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superfluid flow. 

Superfluidity is a very general quantum property of matter at low temperatures, with 
variety of mechanisms and possible nondissipative superfluid currents. These include su- 
percurrent of electric charge in superconductors and mass supercurrent in superfluid ^He, 
where the mechanism of superfluidity is the Cooper pairing; hypercharge supercurrent in the 
vacuum of Standard Model of elementary particle physics, which comes from the Higgs mech- 
anism; supercurrent of color charge in a dense quark matter in quantum chromo-dynamics; 
etc. All these supercurrents have the same origin: the spontaneous breaking of the U{1) or 
higher symmetry related to the conservation of the corresponding charge or particle num- 
ber, which leads to the so called off-diagonal long-range order. This spin supercurrent - 
the superfluid current of spins - is one more representative of superfluid currents. Here 
the U{1) symmetry is the approximate symmetry of spin rotation, which is related to the 
quasi-conservation of spin. 

It appears that the finite life-time of magnons, and non-conservation of spin due to 
the spin-orbital coupling do not prevent the coherence and superfluidity of magnon BEC. 
The non-conservation leads to a decrease of the number of magnons in the Bose gas until 
it disappears completely, but during this relaxation, the coherence of BEC is preserved 
with all the signatures of spin superfluidity: (i) spin supercurrent, which transports the 
magnetization on a macroscopic distance more than 1 cm long; (ii) spin current Josephson 
effect which shows interference between two condensates; (iii) phase-slip processes at the 
critical current; (iv) spin current vortex - a topological defect which is an analog of a 
quantized vortex in superfluids, of an Abrikosov vortex in superconductors, and cosmic 
strings in relativistic theories; (v) Goldstone modes related to the broken f/(l) symmetry - 
phonons in the spin-superfluid magnon gas; etc. 



C. Magnon BEC vs equilibrium magnets 

The magnetic U{1) symmetry is spontaneously broken also in some static magnetic sys- 
tems. Sometimes this phenomenon is also described in terms of BEC of magnons 13Ml6 . 
Let us stress from the beginning that there is the principal difference between the magnetic 
ordering in equilibrium and the BEC of quasiparticles which we are discussing in this review. 

In some magnetic systems, the symmetry breaking phase transition starts when the sys- 
tem becomes softly unstable towards growth of one of the magnon modes. The condensation 
of this mode can be used for the description of the soft mechanism of formation of ferromag- 
netic and antiferromagnetic states (see e.g. However the condensation of this mode 
leads finally to the formation of the true equilibrium ordered state. In the same manner, 
the Bose condensation of phonon modes may serve as a soft mechanism of formation of the 
equilibrium solid crystals [l8[. But this does not mean that the final crystal state is the 
Bose condensate of phonons. 

On the contrary, BEC of quasiparticles is in principle a non-equilibrium phenomenon, 
since quasiparticles (magnons) have a finite life-time. In our case magnons live long enough to 
form a state very close to thermodynamic equilbrium BEC, but still it is not an equilibrium. 
In the final equilibrium state at T = all the magnons will die out. In this respect, the 
growth of a single mode in the non-linear process after a hydrodynamic instability (2l| . 
which has been discussed in terms of the Bose condensation of the classical sound or surface 



waves 



22| . is more close to magnon BEC than equilibrium magnets. 



The other difference is that the ordered magnetic states are states with diagonal long- 
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range order. The magnon BEC is a dynamic state characterized by the off-diagonal long- 
range order (see below), which is the main signature of spin superfluidity. 



II. COHERENT PRECESSION AS MAGNON SUPERFLUID 
A. Spin precession 

The magnetic subsystem which we discuss is the precessing magnetization. In a full 
correspondence with atomic systems, the precessing spins can be either in the normal state 
or in the ordered spin-superfiuid state. In the normal state, spins of atoms are precessing with 
the local frequency determined by the local magnetic field and interactions. In the ordered 
state the precession of all spins is coherent: they spontaneously develop the common global 
frequency and the global phase of precession. 




FIG. 1: The stroboscopic record of the induction decay signal on a frequency about 1 MHz. left: 
During the first stage of about 0.002 s the induction signal completely disappears due to dephasing. 
Then, during about 0.02 s, the spin supercurrent redistributes the magnetization and creates the 
phase coherent precession, which is equivalent to the magnon BEC state. Due to small magnetic 
relaxation, the number of magnons slowly decreases but the precession remains coherent, right: 
The initial part of the magnon BEC signal. 

In pulsed NMR experiments the magnetization is created by an applied static magnetic 
field: M = xH, where x is magnetic susceptibility. Then a pulse of the radio-frequency (RF) 
field Hrf -L H defiects the magnetization by an angle /3, and after that the induction signal 
from the free precession is measured. In the state of the disordered precession, spins almost 
immediately loose the information on the original common phase and frequency induced 
by the RF field, and due to this decoherence the measured induction signal is very small. 
In the ordered state, all spins precess coherently, which means that the whole macroscopic 
magnetization of the sample of volume V is precessing 

+ iMy = M±e'^'+''' , M± = xHV sin p . (2.1) 



6 



This coherent precession is manifested as a huge and long-hved induction signal, Fig. [TJ 
It is important that the coherence of precession is spontaneous: the global frequency and 
the global phase of precession are formed by the system itself and do not depend on the 
frequency and phase of the initial RF pulse. 



B. Off-diagonal long-range order 

The superfiuid atomic systems are characterized by the off-diagonal long-range order 
(ODLRO) [23|. In superfiuid ^He and in the coherent atomic systems the operator of anni- 
hilation of atoms with momentum p = has a non-zero vacuum expectation value: 

(ao) = A/'o/'e*'^*+*" , (2.2) 

where Aq is the number of particles in the Bose condensate, which in the limit of weak 
interactions between the atoms coincides at T = with the total number of atoms M . 

Eq. (12.1 p demonstrates that in the coherent precession the ODLRO is manifested by a 
non-zero vacuum expectation value of the operator of creation of spin: 

{S+) =S, + iSy = :^e^-*+- , (2.3) 

where 7 is the gyromagnetic ratio, which relates magnetic moment and spin. This analogy 
suggests that in the coherent spin precession the role of the particle number M is played 
by the projection of the total spin on the direction of magnetic field Sz- The corresponding 
symmetry group U{1) in magnetic systems is the group of the 0(2) rotations about the 
direction of magnetic field. This quantity Sz is conserved in the absence of the spin-orbit 
interactions. 

The spin-orbit interactions transform the spin angular momentum of the magnetic sub- 
system to the orbital angular momentum, which causes the losses of spin Sz during the 
precession. In our system of superfiuid 3He, the spin-orbit coupling is relatively rather 
small, and thus Sz is quasi-conserved. Because of the losses of spin the precession will fi- 
nally decay, but during its long life time the precession remains coherent. Fig. [1] This is 
similar to the non-conservation of the number of atoms in the laser traps, where the number 
of atoms decreases with time but this does not destroy the coherence of the atomic BEG. 



C. ODLRO and magnon BEC 

The ODLRO in (12. 3 p can be represented in terms of magnon condensation. To view that 
let us use the Holstein-Primakoff transformation, which relates the spin operators with the 
operators of creation and annihilation of magnons 



s — s 

U = aldo = — ^ . (2.6) 
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Eq. (12. 6p relates the number of magnons M to the deviation of spin Sz from its equihbrium 
value ^^(''^"1^1^'''"™) — g — ■y^HV/'^. In the full thermodynamic equilibrium, magnons are 
absent. Each magnon has spin —h, and thus the total spin projection after pumping of M 
magnons into the system by the RF pulse is reduced by the number of magnons, Sz = S—hN". 
The ODLRO in magnon BEC is given by Eq. (12. 2p . where Ao = A/" is the total number of 
magnons (12.61) in the BEC: 

(ao) = ArV2e-*+- = sin ^ e'"*+'" . (2.7) 

Comparing (12. 7p and (12. 2p . one can see that the role of the chemical potential in atomic 
systems /i is played by the global frequency of the coherent precession u, i.e. /i = w. This 
demonstrates that this analogy with the phenomenon of BEC in atomic gases takes place 
only for the dynamic states of a magnetic subsystem- the states of precession. The ordered 
magnetic systems discussed in Refs. are static, and for them the chemical potential 



of magnons is always zero. 

There are two approaches to study the thermodynamics of atomic systems: at fixed parti- 
cle number N or at fixed chemical potential /i. For the magnon BEC, these two approaches 
correspond to two different experimental arrangement: the pulsed NMR and continuous 
NMR, respectively. In the case of free precession after the pulse, the number of magnons 
pumped into the system is conserved (if one neglects the losses of spin). This corresponds 
to the situation with the fixed A^, in which the system itself will choose the global frequency 
of the coherent precession (the magnon chemical potential). The opposite case is the con- 
tinuous NMR, when a small RF field is continuously applied to compensate the losses. In 
this case the frequency of precession is fixed by the frequency of the RF field, n = u = cjrf, 
and now the number of magnons will be adjusted to this frequency to match the resonance 
condition. 

Finally let us mention that in the approach in which M is strictly conserved and has 
quantized integer values, the quantity < ao >= in Eg. (12. 2p . In the same way the quantity 

{Sj^ = in Eq. (12. 3p . if spin Sz is strictly conserved and takes quantized values. This 

means that formally there is no precession if the system is in the quantum state with fixed 
spin quantum number Sz- However, this does not lead to any paradox in the thermodynamic 
limit: in the limit of infinite TV and Sz, the description in terms of the fixed Af (or Sz) is 
equivalent to the description in terms of with the fixed chemical potential fi (or frequency 
u). 



III. PHENOMENOLOGY OF MAGNON SUPERFLUID 
A. Magnon spectrum and magnon mass 

Let us neglect for a moment the anisotropy of spin wave velocity c and the spin-orbit 
interaction. Then the magnon spectrum in ^He-B has the following form: 

where = 'jH. At large momentum, ck ^ u, this spectrum transforms to the linear 
spectrum u = ck of spin waves propagating with velocity c which is on the order of the 
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Fermi velocity vp- At small k, ck <^ uj, this is the spectrum of massive particle 

hk"^ 

Ek = hw{k) , uj{k) =ujl + 7— , , (3.2) 

2m 

where the magnon mass is: 

m = ^ . (3.3) 

Since c ~ fi?, the relative magnitude of the magnon mass compared to the bare mass ma of 
the ^He atom is 

(3.4) 



ms Ep 

where the Fermi energy Ep ~ rrisVp ~ p^p/m^. With the magnon gap hwL ~ 50 at 
ujl ~ IMHz, and Ep ~ IK, one has m ~ lO'^ms. Small mass of these bosons favors the 
Bose condensation. The opposite factor is the small density n of the bosons. From (12. 6p it 
follows that the magnon density is 

n = = ^(l-cos/3) . (3.5) 

n try 

In the typical precessing state of ^He-B the magnon density n is by the same factor hwi/Ep 
smaller than the density of '^He atoms = p^/Sn'^h^ in the liquid 

— . 3.6 

ris Ep 

Nevertheless, one can easily estimate that at temperatures, at which superluid ^He-B exists, 
T ~ 10~^Ep, practically all the magnons must be condensed in the ground state with k = 0, 
i.e. the Bose condensation of magnons must be complete. 

In '^He-B, the typical temperature is big compared to the magnon gap, T ^ hwi- This 
means that according to (13.11) . thermal magnons are mostly the spin waves with linear 
spectrum u{k) = ck and with characteristic momenta kp ~ T/hc. The density of thermal 
magnons is rip ~ kp. At ^He-B temperatures, this density is much smaller then the density 
of condensed magnons and can be neglected. 

The smallness oiup <^T modifies the estimate the temperature of the Bose condensation, 
compared to the atomic gases. Before we start pumping magnons, we have an equilibrium 
system of thermal magnons with /i = 0. After pumping of extra magnons with density n we 
obtain the quasi-equilibrium state in which the number of magnons is temporarily conserved 
and thus the magnon system acquires a non-zero chemical potential, /i 7^ 0. The number of 
extra magnons which can be absorbed by thermal distribution without formation of BEC is 
thus the difference of the distribution function at /x = and /i 7^ at the same temperature: 

n = Y,U{Eu)-f{Ek-li)) . (3.7) 

k 

This quantity reaches its maximum value when ^ = lol- Since ujl '^T one has: 

^-- = ^^2^^-^- (3-8) 

k 
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This gives the dependence of BEC transition temperature on the number of pumped magnons 

At T < Tbec the magnon BEC with A; = must be formed. In ^He-B, Tbec ~ Ep which is 
by 5 to 6 orders of magnitude higher than the temperature at which superfluid ^He exists. 

The above estimate also demonstrates that in sohd state magnetic systems the BEC may 
occur even at room temperature, see Ref. [gI, 0]. 



B. Order parameter and Gross-Pitaevskii equation 

As in the case of the atomic Bose condensates the main physics of the magnon BEC 
can be found from the consideration of the Gross-Pitaevskii equation for the complex order 
parameter. The local order parameter is obtained by extension of Eq. (12.71) to the inho- 
mogeneous case and is determined as the vacuum expectation value of the magnon field 
operator: 

*(r,t) = (^(r,t)\ , n=|^p , Af = I d^r . (3.10) 



where n is magnon density. To avoid the confusion, let us mention that this order parameter 
(I3.10p describes the coherent precession in any system, superfluid or non-superfiuid. It has 
nothing to do with the multi-component order parameter which describes the underlying 
systems - superfluid phases of ^He [24 . 

If the dissipation and pumping of magnons are ignored (on relaxation and pumping 
terms in magnon BEC see Ref. [25j), the corresponding Gross-Pitaevskii equation has the 
conventional form {h= 1): 

dm 5T , , 

-i = , 3.11 
dt 5^1* ^ ' 

where J^{m} is the free energy functional. In the coherent precession, the global frequency 
is constant in space and time (if dissipation is neglected) 

^(r,t) = ^(r)e^'^* , (3.12) 

and the Gross-Pitaevskii equation transforms into the Ginzburg-Landau equation with u = 
/i: 



/^^ = 0. (3.13) 

The important feature of the magnon systems is that that their number density is limited 
2S 

n < Timax = • (3.14) 
For small n ^ ^max the Ginzburg-Landau free energy functional has the conventional form 

-F-/iA/'= j dh + {ul{v) - uj)\^>\^ + F,A . (3.15) 
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Here ^^(r) = 7-ff(r) is the local Larmor frequency, which plays the role of external potential 
U{r) in atomic condensates. The last term Fso(l^'P) contains nonlinearity which come from 
the spin-orbit interaction. It is analogous to the 4-th order term in the atomic BEC, which 
describes the interaction between the atoms. 



C. Spin-orbit interaction as interaction between magnons 

In the magnetic subsystem of superfluid ^He, the interaction term in the Ginzburg-Landau 
free energy is provided by the spin-orbit interaction - interaction between the spin and orbital 
degrees of freedom. Though the structure of superfluid phases of '^He is rather complicated 



and is described by the multi-component superfluid order parameter |2J]), the only output 
needed for investigation of the coherent precession is the structure of the spin-orbit interac- 
tion term Fgod^E'P), which appears to be rather simple. The spin-orbit interaction provides 
the effective interaction between magnons, which can be attractive and repulsive, depending 
on the orientation of spin and orbital orbital degrees with respect to each other and with 
respect to magnetic field. The orbital degrees of freedom in the superfluid phases of ^He 
are characterized by the direction of the orbital momentum of the Cooper pair 1, which 
also marks the axis of the spatial anisotropy of these superfluid liquids. By changing the 
orientation of 1 with respect to magnetic field one is able to regulate the interaction term in 
experiments. 

In superfluid ^He-B, the spin-orbit interaction has a very peculiar properties. The micro- 
scopic derivation (see fl9.13p ) leads to the following form [26]: 

i^so(^, ^,7) = ^^^llisl -l + l cos7(l + s){l + l)r + 
157^ 2 2 

1(1 _ sfil - If + (1 - s^){l - + COS7)] . (3.16) 

It is obtained by averaging of the spin-orbit energy over the fast precession of spins. Here 
s = cos /9, while / = 1 ■ H describes the orientation of the unit vector 1 with respect to the 
direction H of magnetic field. The parameter Ql is the so-called Leggett frequency, which 
characterizes the magnitude of the spin-orbit interaction and thus the shift of the resonance 
frequency from the Larmor value caused by spin-orbit interaction. In typical experimental 
situations, <^ cu^, which means that the frequency shift is relatively small. Finally 7 
is another angle, which characterizes the mutual orientation of spin and orbital degrees of 
freedom. It is a passive quantity: it takes the value corresponding to the minimum of Fgo 
for given s and Z, i.e. 7 = 7(5, 1). 

To obtain Fgod^E'P) in fl3.15p at fixed 1, one must express s via I^E'p: 

I- s = l-cos/3 = ^-^ , (3.17) 

where S = xH/j is spin density. Since Eq. fl3.16p is quadratic in s, the spin-orbit interaction 
contains quadratic and quartic terms in |\E'|. While the quadratic term modifies the potential 
U in the Ginzburg-Landau free energy, the quartic term simulates the interaction between 
magnons. 

The profile of the spin-orbit interaction Fso{s, 1, 7(s, I)) shown in Fig. [2] determines differ- 
ent states of coherent precession and thus different types of magnon BEC in ^He-B, which 
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HPD - magnon BEC 

FIG. 2: The profile of the spin-orbit energy as a function of s = cos /3 and orbital variable /, where 
(3 is the tipping angle of precession and / is the projection of the orbital angular momentum of a 
Cooper pair on the direction of magnetic field. Spontaneous phase-coherent precession emerging 
at / = 1 and s ~ —1/4 is called HPD (Homogeneously Precessing Domain). 

depend on the orientation of the orbital vector 1. The most important of them, which has 
got the name Homogeneously Precessing Domain (HPD), has been discovered 20 years ago 

li,S. 



IV. MAGNON SUPERFLUIDS IN ^HE-B 

A. HPD as unconventional magnon superfluid 

In the right corner of Fig. |2l the minimum of the free energy occurs for / = 1, i.e. 
for the orbital vector 1 oriented along the magnetic field. This means that if there is no 
other orientational effect on the orbital vector 1, the spin-orbit interaction orients it along 
the magnetic field, and one automatically obtains I = cos Pl = 1. The most surprising 
property emerging at such orientation is the existence of the completely flat region in Fig. 
121 The spin-orbit interaction is identically zero in the large range of the tipping angle (3 of 
precession, for 1 > s = cos/3 > —\- 

Fso{(3)i=i = , cos;g>-^, (4.1) 

/ 1 \ 2 1 



FUP)i=i = -^^i{cosP + -] , cos/3<--. (4.2) 
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Using Eq. fl3.17p one obtains the Ginzburg-Landau potential in (13.150 with: 



, \^\^<n, = ^S, (4.3) 

, \^\''>n,= -S. (4.4) 

Eqs. (14. 3 p and (14. 4p demonstrate that when the orbital momentum is oriented along the 
magnetic field, magnons are non-interacting for all densities n below the threshold value 
Uc = {5/4:)S/h. This is a really unconventional gas. 



1^1) 



Fsom) 
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FIG. 3: F — fin for different values of the chemical potential /x = in magnon BEG in ^He-B. 
For fi < U, i.e. for u < ujl, the minimum of -F — /in corresponds to zero number of magnons, 
n = 0. It is the static state without precession. For /i = [/, i.e. for oj = ujl, the energy is the same 
for all densities in the range < n < ric- For fi > U, i.e. for cj > ujl, the minimum of F — fj,n 
corresponds the magnon BEG with density n> Uc- This corresponds to the coherent precession of 
magnetization with tipping angle /3 > 104°. 

The energy profile of F — fin is shown in Fig. |3] for different values of the chemical 
potential fi = u. For /x below the external potential U, i.e. for u < ul, the minimum 
of F — iin corresponds to zero number of magnons, n = 0. It is the static state of ^He-B 
without precession. For fi > U, i.e. for u > ul, the minimum of F — fin corresponds to the 
finite value of the magnon density: 

This shows that the formation of HPD starts with the discontinuous jump to the finite 
density ric = 5S/4h, which corresponds to coherent precession with the large tipping angle 
- the so-called magic Leggett angle, /3c ~ 104° (cos /3c = —1/4). This is distinct from the 
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standard Ginzburg-Landau energy functional (see Eq. (15. 4p for magnon BEG in ^He-A- 
phase below), where the Bose condensate density smoothly starts growing from zero and is 
proportional to n — U for /i > [/. 



B. Two-domain precession 



pumping N magnons magnon BEC formed 



with the same N 



coherent decay of BEC: 

N decreases 
BEC volume decreases 



> 
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magnon 
BEC ^ 
^ domain i'^ ' 
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FIG. 4: Two domains in the coherent precession in ^He-B. lejt. Incoherent spin precession after 
the pulse of the RF held deflects magnetization from its equilibrium value. The total number 
of magnons pumped into the system is A/" = [S — 8^)1^1 magnons. middle: Formation of two 
domains. All M magnons are concentrated in the lower part of the cell, forming the BEC state 
there. The volume of this domain is determined by the magnon density in BEC, V = N /n, where 
n ~ ric- This volume determines the position zq of the domain boundary zq = V/A, where A is the 
area of the cross-section of the cylindrical cell. The position of the interface in turn determines the 
global frequency of presession, which is equal to the local Larmor frequency at the phase boundary, 
fj, = u) = u;l(zq). right Decay of magnon BEC. The number of magnons decreases due to spin and 
energy losses. Since the magnon density in BEC is fixed (it is always close to ric), the relaxation 
leads to the decrease of the volume of the BEC domain. However within this domain the precession 
remains fully coherent. While the phase boundary slowly moves down the frequency of the global 
precession gradually decreases. Fig. [S] left. 



The two states with zero and finite density of magnons, resemble the low-density gas state 
and the high density liquid state of water, respectively. Gas and liquid can be separated in 
the gravitational field: the heavier liquid state will be concentrated in the lower part of the 
vessel. For the magnon BEG, the role of the gravitational field is played by the gradient of 
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magnetic field: 

Vf/ = ViOL = iVH . (4.6) 

Thus applying the gradient of magnetic field along the axis one enforces phase separation, 
Fig. m The static thermodynamic equilibrium state with no magnons is concentrated in 
the region of higher field, where ujl{z) > u, i.e. U{z) > fi. The magnon superfluid - the 
coherently precessing state - occupies the low- field region, where ujl{z) < co, i.e. U{z) < fi. 
This is the so-called Homogeneously Precessing Domain (HPD), in which all spins precess 
with the same frequency u and the same phase a. In typical experiments the gradient is 
small, and magnon density is close to the threshold value ric- 

The interface between the two domains is situated at the position zo where ul{zo) = uj, 
i.e. U (zo) = /i. In the continuous NMR, the chemical potential is fixed by the frequency of 
the RF field: fi = uj-rf, this determines the position of the interface in the experimental cel. 

In the pulsed NMR, the two-domain structure spontaneously emerges after the magneti- 
zation is deflected by the RF pulse (Fig. HJ left and middle). The position of the interface 
between the domains is determined by the number of magnons pumped into the system: 
Af = {S — Sz)/h. The number of magnons is quasi-conserved, i.e. it is well conserved during 
the time of the formation of the two-domain state of precession. That is why the volume of 
the domain occupied by the magnon BEC after its formation is = M jn^. This determined 
the position of the interface, and the chemical potential [i will be agjusted to this position: 

In the absence of the RF field, i.e. without continuous pumping of magnons, the magnon 
BEC decays due to losses of spin. But the precessing domain (HPD) remains in the fully co- 
herent Bose condensate state, while the volume of the magnon superfluid gradually decreases 
due to losses and the domain boundary slowly moves down (Fig. H] right). The frequency 00 
of spontaneous coherence as well as the phase of precession remain homogeneous across the 
whole Bose condensate domain, but the magnitude of the frequency changes with time, since 
it is determined by the Larmor frequency at the position of the interface, fi(t) = UL{z:o{t)). 
The change of frequency during the decay is shown in Fig. [5] left. This frequency change 
during the relaxation was the main observational fact that led Fomin to construct the theory 
of the two-domain precession [9|. 

The details of formation of the magnon BEC are shown in Fig. [H where the stroboscopic 
record of induction decay signal is shown. During the first stage of about 0.002 s the 
induction signal completely disappears due to dephasing. Then, during about 0.02 s, the 
phase coherent precession spontaneously emerges, which is equivalent to the magnon BEC 
state. Due to a weak magnetic relaxation, the number of magnons slowly decreases but the 
precession remains coherent during the whole process of relaxation. The time of formation 
of magnon BEC is essentially shorter than the relaxation time, as clearly shown in Fig. [1] 
right. 

C. London limit: spin and mass supercurrents 

In HPD state of magnon BEC, the magnon density is comparable with the limiting value, 
n ^ (5/8)ninax, as a result the kinetic term in Ginzburg-Landau energy becomes complicated. 
Situation becomes simple in the London limit, where the magnon superfluidity is described 
by the hydrodynamic energy functional written in terms of density n and superfluid velocity 
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Time dependence of induction signal 



lUagnon distribution 




FIG. 5: The amplitude and frequency of the induction decay signal from magnon BEG. left: 
The condensate occupies the domain where the chemical potential ^ > U and radiates the signal 
corresponding to the Larmor frequency at the domain boundary of condensate. With relaxation 
the number of magnons decreases, and the chemical potential moves to the region with a lower 
Larmor frequency, right: The spectroscopic distribution of magnons. Immediately after the RF 
pulse each spin precesses with the local Larmore frequency. After the BEG formation, all the 
spins precess with the common frequency uj and spontaneously emergent common phase a. Due 
to relaxation the number of magnons decreases, leading to the continuously decreasing frequency. 
The small broadening of BEG state is due to relaxation. By comparing the initial broadening of 
the NMR line of about 600 Hz and final broadening of about 0.5 Hz we can estimate that about 
99.9 % of the pumped magnons are in the condensate. 



Superfluidity is phenomenon arising due to spontaneous breaking of U{1) symmetry, 
which in our case is represented by the symmetry group SO (2) of spin rotations about 
direction of magnetic field. In atomic BEC and in helium superfiuids uch symmetry breaking 
leads to a non-zero value of the superfiuid rigidity - the superfiuid density ps which enters 
the non-dissipative supercurrent of particles and thus to mass supercurrent. The same takes 
place for magnon BEC. But since magnons has both mass m and spin —h, they carry both 
the mass and spin supercurrents. The corresponding Goldstone phonon mode of the magnon 
BEC has been experimentally observed: it is manifested as twist oscillations of the precessing 
domain in '^He-B jl^] (see the Section below). 

For n <^ rimax the mass current of magnons is given by the traditional equation: 

J = psVs , Vs = — Va , psiT = 0) = nm . (4.7) 

m 

As in conventional superfiuids, the superfiuid density of the magnon liquid is determined by 
the magnon density n and magnon mass m. To avoid the confusion let us mention that this 
superfiuid density describes the coherent precession in magnetic subsystem and has nothing 
to do with the superfiuid density of the underlying system - the superfiuid ^He. 

Since each magnon carries spin —h, the magnon mass supercurrent is necessarily accom- 
panied by the magnetization supercurrent - the supercurrent of the 2;-component of spin: 
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The nonzero superfluid density > is the main condition for superfluidity. For the HPD 
state of magnon BEC in ^He-B this condition is fulflUed. 

In the London hmit the magnon superfuidity is described by the hydrodynamic energy 
functional. It s similar to that in superfluid liquids and atomic BEC. However, there are im- 
portant differences. One of the is the anisotropy of magnon mass, which leads to anisotropy 
of superfluid density even at T = 0. Another one is the presence of the symmetry breaking 
term which depends explicitly on a. The hydrodynamic energy functional is expressed in 
terms of the canonically conjugated variables - number density n (magnon density) and the 
superfluid velocity Vg, which is expressed via the gradient of the phase of the Bose conden- 
sate a (the phase of the coherent precession of magnetization) . It has the following general 
form: 

P = ^Psij(^)^^si^^sj e(n) - \in\ -Fsb(a, n) . (4.9) 

Here // as before is the chemical potential; and mij{n) is the matrix of magnon masses, which 
in general depends on magnon density n. In ^He-B when the density of magnon increases the 

interaction modifies the magnon mass. It becomes anisotropic. For the magnons propagating 
along the field and in the transverse directions, their mass depends on the tilting angle in 
the following way: 

1 Cm cos/3 (1 — COS/9) 

^T=2^ (4-10) 

mii(n) nwL 

1 cf,(l + cos 3) + c^, (1 - cos B) 

" (4.11) 



m-'-(n) huL 

where the parameters cy and c±_ are on the order of the Fermi velocity vp. 

It is important that anisotropy of magnon mass leads to anisotropic superfluid density 
and superfluid velocity of magnon superfluid: 

Psij{n) = nniij , v^i^h {m~^).. V . (4.12) 

The mass supercurrent becomes isotropic, when it expressed via a: 
dF 

Ji = -— = hnWia , (4.13) 

dvsi 

while the spin supercurrent becomes isotropic, when it is expressed in terms of superfluid 
velocity: 

J^^h{m-')..Jj^hnv,i . (4.14) 

Finally Fgb is the symmetry breaking term which depends explicitly on a. It arises only 
in continuous wave NMR, when the relaxation of magnon BEC is compensated by RF field. 
It gives rise to the mass of Goldstone boson which we discuss later on. 
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D. Supercurrent and linear momentum 



The magnon superfluidity is very similar to superfluidity of the Ai phase of ^He where only 
one spin component is superfluid |2J] : as a result the superfluid mass current is accompanied 
by the superfluid spin current. 

The spin current in magnon superfluids can be obtained directly from the definition of 
linear momentum density in spin systems: 

P = {S - S,)Va = nhVa . (4.15) 

Here we used the fact that S — Sz and a are canonically conjugated variables. This superfluid 
mass current is accompanied by the superfluid current of spins transferred by the magnon 
condensate. It is determined by the spin to mass ratio for the magnon, and because the 
magnon mass is anisotropic, the spin current transferred by the coherent spin precession is 
anisotropic too: 

J: = ^^V.a , (4.16) 

mil [nj 

3] = r^nV±a . (4.17) 

m-^(n) 

The anisotropy of the current in Eqs. f l4.16tH?T7|) is an important modification of the con- 
ventional Bose condensation. This effect is absent in the atomic Bose condensates. 

Let us note that the density of linear momentum of the spin subsystem is not well defined 
globally. While the total momentum of the system is conserved, the canonical momenta of 



the spin and orbital subsystems are not conserved separately [28l-l30|. For the particular 
choice of the linear momentum density in f l4.15p . P is not defined at the points where /3 = vr, 
because at /3 = vr spins are stationary and thus the spin precession angle a is ill defined. This 
is another interesting feature of the magnon superfluids, which becomes important for the 
magnon BEC emerging in normal (non-superfluid) ^He. The latter represents a coherently 
precessing structure at the interface between the equilibrium domain with /3 = and the 
domain with the reversed magnetization, i.e. with /3 = tt [31| . 



E. Goldstone mode of coherent precession — sound in magnon BEC 

In atomic superfluids, sound is the Goldstone mode of the spontaneously broken U{1) 
symmetry. The same sound mode exists in magnon BEC. The sound in magnon subsystem 
has been calculated in Ref.j32| and identified experimentally in Ref. [22]. 

There are several conditions required for the existence and stability of the magnon BEC. 
Some of them are the same as for the conventional atomic BEC, but there are also important 
differences, which we discuss later. One of the conditions is that the compressibility (3m of 
the magnon gas must be positive: 

This condition means that the fourth order term in the Ginzburg-Landau free energy should 
be positive, i.e. the interaction between magnons should be repulsive. The magnon interac- 
tion energy e(n) is provided by spin-orbit (dipole-dipole) interaction. It has a very peculiar 
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form for HPD in ^He-B: 



where 6(x) is Heaviside step function; Ql is Leggett frequency (we assume that Ql ^ i^l)- 
This means that in this state of magnon BEC a stable coherent precession occurs only 
at large enough magnon density n > 5S/Ah, where (PE^o/drt^ > 0. This corresponds to 
cos /3 < —1/4. The magnon BEC state in f l4.3Up also satisfy the condition fl4.18p . while the 
magnon condensates in f l4.29p and in bulk ^He-A are unstable. 

The compressibility of the magnon gas determines the speed of sound propagating in the 
magnon gas. Since the magnons mass is anisotropic the phonon spectrum is also anisotropic: 

■ ■ W P ,^2 771 

(cj ={m ) -^ = n^{m ) . (4.20) 

In the typical experiments with HPD, cos f3 is close to —1/4, i.e. cos /3 = — 1/4 — 0. For such 
/3 one has: 



n d'E.^ 2 f| ^^^^ _ ^ ^^ 21) 



my dn'^ 3 oo^ 



m± dn^ 3 u 
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F. Mass of phonons in magnon superfluid 

As distinct from the conventional superfluids, in magnon BEC one may introduce ex- 
perimentally the symmetry breaking field which smoothly violates the U{1) symmetry and 
induces a small gap (mass) in the phonon spectrum. This mass has been measured. 

The symmetry-breaking term appears in continuous wave NMR, when the relaxation of 
magnon BEC is compensated by RF field. It describes the interaction Fsb(a, n) = — 7Hrf'S 
of the precessing magnetization with the RF field Hrf, which is transverse to the applied 
constant field H. In continuous wave NMR experiments the RF field prescribes the frequency 
of precession, u = wrf, and thus fixes the chemical potential /i; while in the state of free 
precession the chemical potential fi is determined by the number of pumped magnons. The 
symmetry-breaking term depends explicitly on the phase of precession a with respect to the 
direction of the RF- field in the precessing frame: 

F,b = -7i^RF^±cos« = -7i7RP^sin/3(^l-y^ . (4.23) 

Due to explicit dependence on a, this term generates the mass of the Goldstone boson 
(phonon) j33|. For cos/9 = — 1/4 the phonon spectrum becomes: 

^'(k) = {ciy hk, +ml , ml = -^iH^F^ . (4.24) 

V 15 



Two experiments with HPD 34, 35l| reported the gap in the spectrum of the collective mode 

1/2 

of the coherent precession. The measured gap is proportional to H-^p in agreement with 
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G. Spin-current Josephson effect 

The phase coherent precession of magnetizaton in superfluid ^He has all the properties of 
the coherent Bose condensate of magnons. The main spin-superfluid properties of HPD have 
been verified already in the early experiments 20 years ago. These include spin supercurrent 
which transports the magnetization (analog of the mass current in conventional superfluids) ; 
spin current Josephson effect and phase-slip processes at the critical current [36[ [13] (Fig. 
CD. 



H. Spin vortex topological defect of magnon BEC 

Later on the spin current vortex has been observed [SS] - a topological defect which is an 
analog of a quantized vortex in superfluids and of an Abrikosov vortex in superconductors. 
The condensate phase a has 27r winding around the vortex core (see Fig. [S]). The mass 
current around the vortex core is accompanied by spin current. 

In HPD with a single vortex in the central part of cylindrical cell the phase a of precession 
changes by lix around the center, i.e. a is opposite on the opposite sides of the cell. In the 
central part of the cell, i.e. in the vortex core, the magnetization remains vertical and does 
not precess. The magnon BEC with a spin vortex is created by applying the quadrupole 
RF field. For this purpose two parts of the saddle NMR coil are connected in opposite 
directions, so that the phase of RF field (and consequently the phase a) was opposite at the 
opposite sides of the cell. By these NMR coils practically the same HPD signal was observed 
as in the conventional arrangement with the parallel connection of the coils, though with 
a slightly reduced amplitude. This shows that HPD is created with opposite a on opposite 
sides of the cell. To verify this a pair of small pick-up coils are installed at the top of the 
cell connected in usual way. When the RF field is switched off, the pick up coils received 
a very small RF signal from HPD, while the frequency of this signal corresponded to the 
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FIG. 7: Josephson effect for magnon BEC in '^He-B which demonstrates the interference between 
two condensates, left: Spin current as a function of the phase difference across the junction, a2 — ai, 
where ai and 02 are phases of precession in two coherently precessing domains connected by a 
small orifice [right). Different experimental records correspond to a different ratio between the 
diameter of the orifice and the stiffness (magnetic coherence length) of the HPD state. 



full HPD signal. This means that HPD generated the signal with the opposite phase at 
the two sides of the pickup coils, which nearly compensated each other (see Fig. [H]). This 
corresponds to HPD with a circular gradient of a, as shown in Fig. [HI The magnetization is 
oriented vertically in the vortex core. On the periphery of the cell it precesses with tipping 
angle 104°, and with 2tt phase winding around the center. This type of HPD should radiate 
at frequency, which corresponds to the Larmor field on the boundary of HPD, but should 
not produce any signal in the pick-up coil. A small signal appears due to asymmetry of the 
pick-up coil; oscillations of this signal may correspond to nutations of the vortex core. 



I. Critical velocities, coherence length and the vortex core radius 

The phonon spectrum in magnon gas determines the Landau critical velocity of the 
counterflow at which phonons are created. For the case of isotropic sound one has (the 
Landau criterion for the anisotropic case has been derived in [39|): 



vi. = Cs. (4.25) 

For conventional superfluids this suggests that the coherence length and the size of the 
vortex core should be on the order of: 

' ' (4.26) 



mCs ill 

However, for magnon superfluid in the HPD state this gives only the lower bound on the 
core size. The core is larger due to specific profile of the Ginzburg-Landau (dipole) energy 
in Eq.( l4.19p which is strictly zero for cos (3 > —1/4. This leads to the special topological 
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FIG. 8: Spin-mass vortex in magnon BEC. The spin current around the vortex core is accompanied 
by mass current. According to Eq. (j4.26p the size of the vortex core diverges when the HPD domain 
boundary is approached where the local Larmor frequency ul 
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properties of coherent precession (see Ref. [40|). As a result the spin vortex created and 
observed in Ref. ISSj has a continuous core with broken symmetry, similar to vortices in 
superfluid ^He-A [4l|. The size of the continuous core is determined by the proper coherence 
length |42!] which can be found from the competition between the first two terms in the 
Ginzburg-Landau free energy in Eq.f l3.15|) : 



h 



This coherence length determines also the critical velocity for creation of vortices: 

h 



(4.27) 



(4.28) 



For large tipping angles of precession the symmetry of the vortex core is restored: the 
vortex becomes singular with the core radius rcore ^ c/Ql j43j . 

Other topological defects possible in the coherent precession beyond the Ginzburg-Landau 
model of magnon BEC are discussed in j4o|. In BEC of excitations, the topological defects 
have been detected in exciton-polariton condensate 44|. Among them the half quantum 
vortices - vortices with half of the circulation quantum. Such vortices are topologically 
stable in the superfluid ^He-A t45|], but they still remain elusive there. 
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FIG. 9: NMR signature of the spin vortex in magnon BEC (HPD). The frequency decays with 
time in the same way as in conventional vortex-free precession, while the amplitude of the HPD 
signal is nearly zero because of compensation of signals from opposite sides of the cells where the 
phase a of precession differs by tt. 

J. Other states of magnon BEC in ^He-B 



Recent experiments in ^He-B allowed to probe the BEC states that emerge in the valley 
on the other side of the energy barrier in Fig. |2l This became possible by immersing the 
superfluid ^He in a very porous material called aerogel. By squeezing or stretching the 
aerogel sample, one creates the global anisotropy which captures the orbital vector 1. This 
allows to orient the orbital vector 1 in the desirable direction with respect to magnetic field 



46, 47 



For the transverse orientation of 1, i.e. for I = 0, two new BEC states have been identified. 
One of them exists at < S/h and has the following form of spin-orbit interaction 
obtained from Eq. f l3.16p (we omit for simplicity the constant term): 



X 



VP' 



1^ 



1^1 



S 



(4.29) 



This state has an attractive interaction between magnons, and is unstable since the com- 
pressibility Pm of the magnon gas in (14.181) is negative: cPe/drP < 0. 
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The other state exists at > S/h and has the following form of spin-orbit interaction 

2 



fin ^ - 2 



> 



2072 ^ V 

This state has repulsive interaction between magnons and is stable. 



(4.30) 



V. MAGNON BEG IN ^HE-A 

As in the case of ^He-B, all the information on the ^He-A order parameter needed to 
study the coherent precession is encoded in the spin-orbit interaction. 



A. Instability of magnon BEG in bulk ^He-A 

For ^He-A, the spin-orbit interaction averaged over the fast precession has the following 
form 48 : 
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(5.1) 



In a static bulk ^He-A, when \& = 0, the spin-orbit energy Fg^ in Eq. fl5.ip is minimized when 
the orbital vector 1 is perpendicular to magnetic field, i.e. for I = 0. Then one has 



Fso(|^M = o) 
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(5.2) 



with a negative quartic term. The attractive interaction between magnons destabilizes 
the BEC, which means that homogeneous precession of magnetization in ^He-A becomes 
unstable. This instability predicted by Fomin jioj was experimentally confirmed 50 . 
However, as follows from (15. ip . at sufficiently large magnon density n = |\E'|2 



8 + \/8 



S > n> 
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(5.3) 



the factor in front of becomes negative. Therefore it becomes energetically favorable to 
orient the orbital momentum I along the magnetic field, / = 1. For this orientation one 
obtains the Ginzburg-Landau free energy with 



i^so(|^M = l) 
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(5.4) 



It corresponds to the conventional Ginzburg-Landau free energy in atomic BEC. The 
quadratic term modifies the potential U; the quartic term is now positive. 

In the language of BEC, this means that, with increasing the density of Bose conden- 
sate, the originally attractive interaction between magnons should spontaneously become 
repulsive when the critical magnon density Uc = S{8 — \/8)/7 is reached. If this happens, 



24 



-0,2 



-0,4 



-0,6 



-0,8 



— , — 1 — 1 — 1 — 1 — 1 — 1 — 1 — I — 

i[i=^;i-u>o 
































|U = ^- 


U = 





0,5 



1.5 

^^/S = 1 - cosp 



FIG. 10: F — fj,n for different values of the chemical potential fi > U in magnon BEC in ^He-A. 
Magnon BEC in ^He-A is similar to BEC in atomic gases. 



the magnon BEC becomes stable and in this way the state with spontaneous coherent pre- 
cession could be formed [48[. This self-stabilization effect is similar to the effect of Q-ball, 
where bosons create the potential well in which they condense (we shall discuss the Q-ball 
phenomenon in magnon BEC later on). However, such a self-sustaining BEC with originally 
attractive boson interaction has not been achieved experimentally in bulk ^He-A, most prob- 
ably because of the large dissipation, due to which the threshold value of the condensate 
density has not been reached. 



B. Magnon BEC of ^He-A in deformed aerogel 

Finally the fixed orientation of the orbital vector 1 has been achieved in ^He-A confined in 
aerogel - the material with high porosity, which is about 98% of volume. Silicon strands of 
aerogel play the role of impurities with local anisotropy along the strands. According to the 
Larkin-Imry-Ma effect, the random anisotropy suppresses the orientational long-range order 
of the orbital vector 1; however, when the aerogel sample is deformed the long-range order 
of i is restored [5l|. Experiments with globally squeezed aerogel 52| demonstrated that a 
uni-axial deformation by about 1% is sufficient for global orientation of the vector 1 along 
the anisotropy axis. When magnetic field is also oriented oriented along the anisotropy axis 
one obtains the required geometry with / = 1, at which the magnon BEC in ^He-A becomes 
stable. The first indication of coherent precession in ^He-A has been reported in 53|, |54 
and confirmed in (Hsf. Contrary to the unconventional magnon BEC in the form of HPD in 
^He-B, the magnon BEC emerging in the superfluid ^He-A is in one-to-one correspondence 
with the atomic BEC, see Fig. (TDl For fi > U, the condensate density determined from 
equation dF/ dn = fi continuously grows from zero as oc /i — f/. 
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For / = 1 the Ginzburg-Landau free energy acquires the standard form: 

F = Jd'r (^-^ + K(r) - /i)|vl/p + ^b\^\'^ , (5.5) 
where we modified the chemical potential by the constant frequency shift: 

= . (5.6) 

and the parameter b of repulsive magnon interaction is 

■ ' (5.7) 



2ujS 

At fi > ojl, magnon BEG must be formed with density 

i^r = ^^- (5.8) 

This is distinct from ^He-B, where condensation starts with finite condensate density. Eq. 
(15. 8 p corresponds to the following dependence of the frequency shift on tipping angle (3 of 
coherence precession: 

CO - uj^ = --J± cos 13 . (5.9) 
The final proof of the coherence of precession in ^He-A in aerogel was the observation 



of the free precession after a pulsed NMR and also after a switch off the GW NMR j55 
In conclusion, in contrast to the homogeneously precessing domain (HPD) in '^He-B, the 
magnon Bose condensation in ^He-A obeys the standard Gross-Pitaevskii equation. In bulk 
^He-A, the Bose condensate of magnons is unstable because of the attractive interaction 
between magnons. In ^He-A confined in aerogel, the repulsive interaction is achieved by the 
proper deformation of the aerogel sample, and the Bose condensate becomes stable. 



VI. MAGNON BEG IN MAGNETIG TRAP AND Q-BALL 

There are many new physical phenomena related to the Bose condensation of magnons, 
which have been observed after the discovery of HPD. These include in particular compact 
objects with finite number of the Bose condensed magnons. At small number N of the 
pumped magnons, the system is similar to the Bose condensate of the ultracold atoms in 
harmonic traps, while at larger N the analog of the Q-ball in particle physics develops j56| . 

A Q-ball is a non-topological soliton solution in field theories containing a complex scalar 
field \E'. Q-balls are stabilized due to the conservation of the global U{1) charge Q (sTj : 
they exist if the energy minimum develops at nonzero at fixed Q. At the quantum level, 
Q-ball is formed due to suitable attractive interaction that binds the quanta of \E'-field into 
a large compact object. In some modern SUSY scenarios Q-balls are considered as a heavy 
particle-like objects, with Q being the baryon and/or lepton number. For many conceivable 
alternatives, Q-balls may contribute significantly to the dark matter and baryon contents of 
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FIG. 11: The trapping potential Eq. (16.2p is formed in the cyhndrically symmetric "flare-out" 
texture of the orbital angular momentum 1 (dashed lines) in a shallow minimum of the vertical 
magnetic field {right). The arrows represent the precessingmagnetization M, which precesses 



coherently within the condensate droplet (dark blue). From [63|] 



the Universe, as described in revie w |58l| . Stable cosmological Q-balls can be searched for in 
existing and planned experiments [59|. 

The Q-ball is a rather general physical object, which in principle can be formed in con- 
densed matter systems. In particular, Q-balls were suggested in the atomic Bose-Einstein 
condensates 60|]. In ^He-B, the Q-balls are formed as special states of phase coherent pre- 
cession of magnetization. The role of the Q-charge is played by the projection of the total 
spin of the system on the axis of magnetic field, which is a rather well conserved quantity at 
low temperature, or which is the same the magnon number A/". At the quantum level, this 
Q-ball is a compact object formed by magnons - quanta of the corresponding \l/-field. 

In ^He-B the Q-balls are formed at low temperatures, when homogeneous magnon BEC in 
the form of Homogeneously Precessing Domain (HPD) becomes unstable due to parametric 



Suhl instability [6l|, |62| , which we shall discuss later. At low temperatures the condensate 
can be formed only in a trap, similar to that in atomic gases ll|], and the Q balls are formed 
in traps. 



A. Magnetic trap in ^He-B 

A cylindricall y sy mmetric magnetic trap for magnon BEC in ^He-B is schematically 
shown in Fig. [11] |63| . The confinement potential U\\ (z) in the axial direction is produced by 
a small pinch coil, while in the radial direction the well U±{r) is formed by the cylindrically 
symmetric flare-out texture of the orbital vector 1: 

2Q2 

f/(r) = U\\{z) + U^{r) = + ^ (1 - ^(0) , (6-1) 

where as before / = 1 ■ H describes the orientation of the unit vector 1 with respect to 
the direction H of magnetic field. On the side wall of the cylindrical container the orbital 
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momentum 1 is normal to the wall, while in the center it is parallel to the axially oriented 
applied magnetic field. This produces a minimum of the potential U±{r) on th e cy linder 
axis. Close to the axis the angle (3l varies linearly with distance r from the axis [4l|. As a 
result, the potential U (r) reduces to that of a usual harmonic trap used for the confinement 
of dilute Bose gases 

U{r) = U{0) + '^{u:y + coy) , (6.2) 

The oscillator frequency Ur of the well in the radial direction can be adjusted by applying 
rotation, since adding vortex- free superfiuid fiow and/or rectilinear vortex lines modifies the 
fiare-out texture. Low-amplitude standing spin waves in the trap have the conventional 
spectrum 

= '^l(O) + Ur{m. + 1) + w^(n + 1/2) . (6.3) 
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FIG. 12: Numerical simulation of the Q-ball growing in the center of the cell. At small Q, when 
the /3-angle of spin S is less than about 10°, the l-texture is practically Q-independent (see /3l(x) 
marked by solid line). At larger Q the potential well for magnons becomes wider: the dashed line 
corresponds to the l-texture formed by Q-ball with (3 = 20° in the middle of the Q-ball; and the 
dotted line - with /3 = 33° . For details see [65 ] . 



B. Ground-state condensate and self-localization 

Typically the condensate is formed in the ground state (0,0) of the trap. When the 
number of magnons A/" in the ground state increases, the interactions between them become 
important and the condensate wave function starts deviating from the Gaussian form of an 
ideal gas. As distinct from a system of cold atoms, the peculiarity of the magnon Ginzburg- 
Landau functional in Eq. (I3.16p is that the prefactor of the quartic term, which describes 
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the interactions between the magnons, is not a constant, but oc w^r^ and thus is small in 
the region of the trap and can be neglected. 

Under conditions of our experiment the main effect the self-localization discussed in f56|. 
At high density of magnons they start to influence the radial 1-texture. According to f l3.16p 
the condensate leads to the preferable orientation of 1 parallel to magnetic field, / = 1, in 
the region of the trap. As a result the potential well becomes wider and the energy of the 
level in the trap decreases. This allows to incorporate more magnons at this same level by 
sweeping the frequency up. This is equivalent to effective attractive interaction induced by 
the exchange of the quanta of the I-field. 

In the language of relativistic quantum fields, the self-localization is caused by interaction 
between the charged and neutral fields, where the neutral field provides the potential for the 
charged one. In the process of self-localization the charged field modifies locally the neutral 
field so that the potential well is formed in which the charge is condensed. Such a trapped 
condensate resembles a droplet of quantum field with a fixed value of field quanta Q - the 
Q-ball jH^]. The charge Q corresponds to the spin or, equivalently, to the magnon number 
M. Experimentally the Q-ball in '^He-B ^] is manifested as l ong -lived ringing (of up to an 
hour!) of the free induction decay after a NMR tipping pulse [GJ]. In steady state it can be 



maintained by CW RF pumping (Fig. [T3ll . and even by off- resonance excitation [65|, |66 . 

One can estimate the radius of the self-localization as a function of the magnon number 
at large U |63|: 



R{N)^ ^-(at) ' -^^-^^ ^6.4) 

where is the harmonic oscillator length in the original (non-modified) radial trap, and Mc 
is the characteristic number at which interaction becomes important. Assuming that the 
trap is flat in the region of the condensate, /3l = 0, one obtains p = 1/2. For comparison, for 
an atomic BEC with repulsive inter-particle interaction the exponent in the Thomas-Fermi 
limit is p = 1/5 [ll|. Compare also the change in chemical potential due to interaction in 
the atomic BEC in the Thomas-Fermi regime, and the frequency shift in magnon BEC: 

/i-f/(0) -^j , atomic BEC, (6.5) 



-1 



uj — a;L(0) ~ uJr j ' "^cLgnon BEC . (6.6) 

With the growing Q-ball, its frequency uj decreases approaching the Larmor frequency 
asymptotically. As distinct from the atomic BEC, the magnon BEC droplet has d^/dM < 0. 
This determines the way in which the magnon condensate is grown in a cw NMR measure- 
ment, as seen in Fig. [13] for the ground state of the trap in Fig. [11] The magnons are created 
when the frequency u of the applied RF field is swept down and crosses the ground state 
level Uqq. When u is reduced further, the number of magnons s follows from Eq. f l6.6p . 

Q-balls can be also formed in pulsed NMR measurements. As distinct from cw NMR, 
where the Q-balls are generated starting continuously from Af = 0, in pulsed NMR a Q- 
ball is formed after a large is pumped to the cell. In this case the 3D Q-ball is often 
formed on the axis of the fiared-out texture, away from the horizontal walls J88|]. This clearly 
demonstrates the effect of self-localization: the main part of the pumped charge relaxes but 
the rest of J\f starts to concentrate at some place on the axis, digging a potential well there 
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FIG. 13: Formation of the magnon condensate droplet in the ground state of the trap in a cw NMR 
measurement. A4± grows while the frequency of applied RF field / = u}/{2tt) is swept down (or 
equivalently, the steady magnetic field H is swept up) , in agreement with Eq. (j6.6p . The condensate 
is destroyed (vertical lines), when energy dissipation exceeds the RF pumping and therefore the 
maximum volume of the condensate droplet grows with increasing RF excitation. 

and attracting the charge from other places of the container. Moreover, this also shows 
that Q-balls are not necessarily formed at the bottom of the original potential: Q-ball may 
dig the potential well in a different place. In the formation of Q-ball with off-resonance 
excitation 



66|, 

the effect of self-localization also plays a crucial role. 
In conclusion, Q-balls represent a new phase coherent state of Larmor precession. They 
emerge at low T, when the homogeneous bulk BEC of magnons (HPD) becomes unstable. 
These Q-balls are compact objects which exist due to the conservation of the global U{1) 
charge Q = Sz- At small Q they are stabilized in the potential well, while at large Q the 
effect of self-localization is observed. In terms of relativistic quantum fields the localization 
is caused by the peculiar interaction between the charged and neutral fields 67| • The neutral 
field / provides the potential for the charged field the charged field modifies locally the 
neutral field so that the potential well is formed in which the charge Q is further condensed. 



C. Non-ground-state condensates 

The condensate can be also formed when one starts filling magnons to one of the levels 
(m,n) in Eq. fl6.3p . Then one obtains the non-ground-state condensate 63|. The formation 
of non-ground-state condensates has been proposed for cold atoms [gsI], but as a dynamic 
mixture of the ground state and an excited level. It was suggested to use resonant modulation 
of either the trap potential or the atomic scattering length, eg. by applying a temporal 
modulation of the atomic interactions via the Feshbach resonance technique. In contrast to 
such schemes with atomic condensates in optical traps, the properties of magnon condensates 
make it possible to populate different excited trap levels by pumping magnons resonantly 
directly to these levels. This fact shows that condensates formed from excitations lead to 
interesting new consequences. 
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FIG. 14: Twisted core of non-axisymmetric vortex in ^He-B. The gradient of the Goldstone field 
V(/) along the string corresponds to the superconducting current along the superconducting cosmic 
string. Such twisted core has been obtained and detected using coherent precession of magnetiza- 



tion 



The excited states (m,n) can be populated [63|, because their frequencies Wmn also de- 
crease with increasing the magnon number Af. The condensate in the state (m,n) grows 
when the frequency of RF field is swept down and crosses the level Wmn(O) from above (63| . 
An excited state is populated only if magnons are pumped directly to this level by apply- 
ing an RF field with the exact frequency Ujani-^)- The condensate in the excited state is 
metastable: it is supported by continuous pumping at WmnlA/") or otherwise decays to the 
ground state cuqo- Thus in pulsed NMR the long-lived free precession, which is still ringing 
well after the RF tipping pulse, is produced by the ground state population. The spectrum 
f l6.3p and its evolution due to self-localization can be conveniently examined experimentally 
by modifying the trap in Fig. [TTl by applying rotation or by changing the depth of the axial 
field minimum. 



VII. EXPLOITING BOSE CONDENSATE OF MAGNONS 

The Bose condensation of magnons in superfluid ^He-B has many practical applications. 
In Helsinki, owing to the extreme sensitivity of the Bose condensate to textural inhomogene- 
ity, the phenomenon of Bose condensation has been applied to studies of supercurrents and 
topological defects in ^He-B in rotating cryostat. The measurement technique was called 



HPD spectroscopy |69|, [70 



A. Observation of Witten string 



In particular, HPD spectroscopy provided direct experimental evidence for broken axial 
symmetry in the core of quantized vortices in ^He-B 71 . 



The dominating area of the phase diagram of the vortex states in ^He-B is occupied by 
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vortices with non-axisymmetric cores, i.e. vortices with the spontaneously broken rotational 



SO (2) of the core calculated in Refs. j72|, |73||. The core with broken rotational symmetry 



can be considered as a pair of half-quantum vortices, connected by a non-topological soliton 
wall (see Fig. [T^ . The separation of the half-quantum vortices increases with decreasing 
pressure and thus the double-core structure is most pronounced at zero pressure. 

In the physics of cosmic strings, an analogous breaking of continuous symmetry in the 
core was first discussed by Witten [t^I, who considered the spontaneous breaking of the 
electromagnetic gauge symmetry U{1). Since the same symmetry group is broken in the 
condensed matter superconductors, one can say that in the core of the cosmic string there 
appears the superconductivity of the electric charges, hence the name 'superconducting 
cosmic strings'. 

For the ^He-B vortices, the spontaneous breaking of the S0{2) symmetry in the core 
leads to the Goldstone bosons - the mode in which the degeneracy parameter, the axis of 
anisotropy of the vortex core, is oscillating. The homogeneous magnon condensate, the HPD 
state, has been used to study the structure and twisting dynamics of the non-axisymmetric 
core of the low-temperature vortex in '^He-B [7l| . This is because the coherent precession of 
magnetization excites the vibrational Goldstone mode via spin-orbit interaction. Moreover, 
due to spin-orbit interaction the precessing magnetization rotates the core around its axis 
with constant angular velocity. In addition, since the core was pinned on the top and the 
bottom of the container, it was possible even to screw the core (see Fig. [H]). Such a twisted 
core corresponds to the Witten superconducting string with the electric supercurrent along 
the core. The rigidity of twisted core differs from that of the straight core. This allowed 
a detailed study of the Goldstone mode of the vortex core resulting from the spontaneous 
violation of rotational U{1) symmetry in the core (7l| . 

B. Observation of spin-mass vortex 

There are different types of the topological defects in the (non-precessing) ^He-B. Among 
them there is a Z2 spin vortex - topological defect of the order parameter matrix Rai in 
f l9.8p . Due to spin-orbit coupling this defect serves as the termination line of the topological 
soliton wall, and due to the soliton tension it cannot be stabilized in the rotating vessel. 
However, spin and mass vortices attract each other and form the combined defect, the so- 
called spin- mass vortex which can be stabilized under rotation (see Fig. [T5|) . The spin-mass 
vortices also form molecules where the soliton serves as chemical bond. These defects - 
spin-mass vortex connected with the wall by soliton and bound pairs of spin-mass vortices 



- have been observed and studied using HPD spectroscopy [75 



C. Magnon condensates in aerogel 

HPD spectroscopy proved to be extremely useful for the investigation of the superfiuid 



order parameter in a novel system - superfiuid He confined in aerogel [53|, I76l-l80 



D. Towards observation of Majorana fermions 

The condensate in the magnetic trap can be used to continue measurements down to 0.1 - 
0.2 Tc. Of great current interest are the fermion states in the vortex core, especially the still 
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FIG. 15: Vortices in rotating ^He-B. Mass vortices form a regular structure like Abrikosov vor- 
tices in an applied magnetic field. If the number of vortices is less than equilibrium number for 
given rotation velocity, vortices are collected in the vortex cluster. Within the cluster the average 
superfluid velocity (vg) = Vn. On the periphery there is a region void of vortices - the counterflow 
region, where (vg) 7^ Vn. Spin-mass vortices can be created and stabilized in the rotating vessel. 
The confining potential produced by the soliton wall is compensated by logarithmic repulsion of 
vortices forming the vortex pair - the doubly quantized vortex inside the cluster. A single spin- 
mass vortex is stabilized at the periphery of the cluster by the combined effect of soliton tension 
and Magnus force acting on the mass part of the vortex. 



elusive Majorana fermions. Recently the topologically nontrivial gapless and gapped phases 
of matter, topological insulator states in semimetals, superconductors, and superfluids have 
attracted attention. ^He-B is the best representative of a 3- dimensional coherent quantum 
system with time reversal symmetry. Its nontrivial topology gives rise to gapless Andreev- 
Majorana fermions as surface states There is now experimental evidence for Andreev 

surface states in ^He-B at a solid wall 0, 85|, but the Majorana signature of these fermions 



- the linear 'relativistic' energy spectrum at low energy - can be observed only at extreme 
low temperatures, when the thermal quasiparticles in bulk ^He-B are exponentially depleted. 
Majorana fermions, both in the vortex core and at surfaces, are then expected to give the 
main contribution to thermodynamics and dissipation, with a power-law dependence of the 
physical quantities on T. In rotation or by moving the magnon condensate droplet next to 
the wall, one will be able to probe Majorana fermion states. 



VIII. MAGNON BEG IN OTHER SYSTEMS 



A. Magnon BEG in normal '^He 



A very long lived induction signal was observed in normal Fermi liquids: in spin-polarized 



■^He-^He solutions 86| and in normal liquid ^He [87|. It was explained as a coherently 
precessing structure at the interface between the equilibrium domain and the domain with 
the reversed magnetization 3l|]. It would be interesting to treat this type of dynamic 
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magnetic ordering as a new mode of magnon BEC. 



B. Magnons condensation in yttrium-iron garnet 

In all the magnon systems under consideration, the typical temperature is big compared 
to the magnon gap, T ^ hwi. In ^He-B this means that according to fl3.ll) thermal magnons 
are spin waves with linear spectrum uj{k) = ck, with characteristic momenta ~ T/hc. 
The density of thermal magnons nx ~ kj^ is much smaller than the density of condensed 
magnons and in ^He-B can be neglected. 

In yttrium-iron garnet (YIG) situation is opposite: the temperature is high and number 
density of the magnons concentrated at small momentum is small compared to the thermal 
magnons, n <C n-r [sl, 0]- Magnons in yttrium-iron garnet (YIG) films have the quasi 2D 
spectrum: 0, 0] 

u^ik,, A:,) = A„ + + I^^LZM! , (8.1) 
Zrriy lirix 

where magnetic field is along x\ the gap in the lowest branch Aq = 2.1 GHz = 101 mK 
at = 700 Oe Q and Aq = 2.9GHz at = 1000 Oe 3; the position of the minimum 
kQ = 5-10^ 1/cm [i]; the anisotropic magnon mass can be probably estimated as m^, ~ /cq/A 
with my being somewhat bigger, both are of order of electron mass. 

In two-dimensional systems the number of extra magnons - the difference of the distri- 
bution function at = and 7^ - is determined by low energy Rayleigh- Jeans part of 
the spectrum: 



n 



If one neglects the contribution of the higher levels and consider the 2D gas, the Eg. (18. 2 p 
becomes 

^extra = -./ma,m„ lu , (8.3) 

In 2D, all extra magnons can be absorbed by thermal distribution at any temperature 
without formation of Bose condensate. The larger is the number n of the pumped magnons 
the closer is /i to Aq, but yU never crosses Aq. At large n the chemical potential exponentially 
approaches Aq from below and the width of the distribution becomes exponentially narrow: 

{5kyf {5kxf Ao - / 2TxhnM 



AAA exp - \ . (8.4) 

If one uses the 2D number density n = 6N d with the film thickness ci = 5 /im and 3D 
number density 5N ~ 5 ■ 10^^ cm~'^, one obtains that at room temperature the exponent is 



271 Hum 2ti /S.q 5N d ^ ^ 



Ty/mxTriy T IO/cq 



10^ . (8.5) 



If this estimation is correct, the peak should be extremely narrow, so that all extra magnons 
are concentrated at the lowest level of the discrete spectrum. However, there are other 
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contributions to the width of the peak due to: finite resolution of spectrometer, magnon 
interaction, finite hfe time of magnons and the infiuence of the higher discrete levels n 7^ 0. 

In any case, the process of the concentration of extra magnons in the states very close to 
the lowest energy is the signature of the BEC of magnons. The main property of the room 
temperature BEC in YIG is that the transition temperature is only slightly higher than 
temperature, Tc — T ^ T; as a result the number of condensed magnons is small compared 
to the number of thermal magnons: n <^ ut- Situation with magnon BEC in ^He is the 
opposite, one has T <^Tc and thus n ^ ut- 



C. Magnon BEC vs planar ferromagnet 

The spin density in coherently precessing state HPD state of ^He-B and in YIG are 
correpondingly 

+ iSy = ^^6^^*+^" , (8.6) 

and 

+ iSy = cos(A;ox)e^'^*+^" . (8.7) 



On the other hand in the equilibrium planar ferromagnet one has |17l . |89 



+ iSy = S^e''' . (8.8) 

This means that in principle, the equilibrium planar ferromagnet can be also described in 
terms of the ODLRO. 

Magnons were originally determined as second quantized modes in the background of 
stationary state with magnetization along z axis. Both the static state in Eq. 08.81) and 
precessing states (18. 6p and (18. 7p can be interpreted as BEC of these original magnons. 
On the other hand, the stationary and precessing states can be presented as new vacuum 
states, the time independent and the time dependent vacua respectively. The excitations - 
phonons - are the second quantized modes in the background of a new vacuum. What is 
the principle difference between the stationary vacuum of planar ferromagnet and the time 
dependent vacuum of coherent precession? 

The major point which distinguishes the HPD state (18.61) in ^He-B and the coherent 
precession (18. 7p in YIG from the equilibrium magnetic states is the conservation (or quasi- 
conservation) of the U{1) charge Q. The charge Q is played by the spin projection Sz in 
magnetic materials (or the related number N of magnons) and by number of atoms in 
atomic BEC. This conservation gives rise to the chemical potential fi = dE/dN, which 
is the precession frequency u in magnetic systems. On the contrary, the static state in 
Eq. (l8.8p does not contain the chemical potential u, i.e. the conservation is not in the origin 
of formation of the static equilibrium state; the chemical potential of magnons in a fully 
equilibrium state is always strictly zero, /i = w = 0. 

The spin-orbit interaction violates the conservation of Sz, as a result the life time of 
magnons is finite. For the precessing states (18.60 and (18. 7p this leads to the finite life time 
of the coherent precession. To support the steady state of precession the pumping of spin 
and energy is required. On the contrary, the spin-orbit interaction does not destroy the 
long-range magnetic order in the static state (18.80 : this is fully equilibrium state which does 
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not decay and thus does not require pumping: the hfe time of static state is macroscopically 
large and thus by many orders of magnitude exceeds the magnon relaxation time. That is 
why a planar ferromagnet flS.Sp is just one more equilibrium state of quantum vacuum, in 
addition to the easy axis ferromagnetic state, rather than the magnon condensate. 

The property of (quasi) conservation of the U{1) charge Q distinguishes the coherent 
precession from the other coherent phenomena, such as optical lasers and standing waves. 
For the real BEC one needs the conservation of particle number or charge Q during the time 
of equilibration. BEC occurs due to the thermodynamics, when the number of particles (or 
charge Q) cannot be accommodated by thermal distribution, and as a result the extra part 
must be accumulated in the lowest energy state. This is the essence of BEC. 

Photons and phonons can also form the true BEC (in thermodynamic sense) under pump- 
ing, again if the lifetime is larger than thermalization time. These BEC states are certainly 
different from such coherent states as optical lasers and from the equilibrium deformations 
of solids. 



IX. BEYOND MAGNON BEC: SUHL INSTABILITY OF COHERENT PRECES- 
SION 

A. Catastrophic relaxation of magnon BEC 

The instability of BEC which we discuss here is applicable only to the BEC of magnon 
quasiparticles, and is irrelevant for atomic condensates. For quasiparticles the U{1) sym- 
metry is not strictly conserved. For magnetic subsystem, this is the 5*0(2) symmetry with 
respect to spin rotations in the plane perpendicular to magnetic field, and it is violated by 
spin-orbit interactions. The magnon BEC is a time dependent process, and it may experience 
instabilities which do not occur in equilibrium condensates of stable particles. In 1989 it was 
found that the original magnon condensate - the HPD state - abruptly looses its stability 
below about 0.4 Tc [61,]. This was called the catastrophic relaxation. This phenomenon was 
left unexplained for a long time until the reason was established: in the low-temperature 
regime, where dissipation becomes sufficiently small, the Suhl instability destroys the homo- 
geneous precession [90l-9^. This is the parametric instability, which leads to decay of HPD 



due to the parametric amplification of spin wave modes. The latter modes are different from 
the magnons which we discuss here: they represent another branch of the collective modes 
of superfluid ^He-B. The instability occurs because the spin-orbit interactions violate the 
U{1) symmetry. 

For magnetically ordered systems the instability of homogeneous precession is a well 



known phenomenon. Suhl [93|] explained it in terms of parametric instability of the mode 
of precession with respect to excitations of pairs of spin waves satisfying the condition of 
resonance: 

nuL = Ws(k) + w^(-k) , (9.1) 



where ool is the precession frequency and n is integer (see also review [9J]). All the mag- 
netic systems, where the Suhl instability has been observed, are anisotropic. In particular, 
as quantum solids and liquids are concerned, Suhl instability has been observed in anti- 



ferromagnetic solid ^He [95|, and has been predicted by Fomin for anisotropic superfluid 
liquid ^He-A 96] where it has been observed later [s^]- Due to the extreme isotropy of 
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^He-B and due to the unique symmetry of the spin-orbit interaction, the Suhl instabihty 
was not expected there. 

However, under conditions of the experiment, the boundary conditions on the wall of 
container induce the texture of the order parameter in which the orbital vector 1 deviates 
from its symmetric orientation along the magnetic field H. The symmetry of the spin-orbit 
interaction is violated providing the additional term in the interaction between the m odes, 
which is dominating in typical experiments with the catastrophic relaxation 98- lOOj . 



B. Precessing states and their symmetry 

To describe the interaction of magnon condensate with the other modes of superfiuid '^He- 
B, we must go beyond the magnon BEC description and consider all the degrees of freedom 
of homogeneous free precession in external magnetic field H. In liquid '^He the spin-orbit 
(dipole-dipole) interaction is weak. If it is neglected, we can apply the powerful Larmor 
theorem, according to which, in the spin-space coordinate frame rotating with the Larmor 
frequency the effect of magnetic field on spins of the ^He atoms is completely compensated. 
This follows from the observation that in dynamics the time derivative of a spin vector enters 
equations together with the Larmor frequency vector cj^ = 7H 

Af = 5<f-a;i X f. (9.2) 

In other words, the Pauli magnetic field acts on spin vectors as time component of the 
effective SO (3) gauge field, 

Ao = 7H. (9.3) 

We shall use this equation later in Sec. |X]for discussion of spin currents and spin-quantum 
Hall effect. 

Since the magnetic field becomes irrelevant, the symmetry group of the physical laws 
in the precessing frame becomes the same as in the absence of the field. If the spin-orbit 
interaction is neglected, it is the product of the SO^ group of orbital rotations and the SO^ 
group of spin rotations: 

G = SO^ X sol , (9.4) 

The only difference is that while 5*03 is the group of orbital rotations in the laboratory 
frame, the SOl rotations are in the precessing frame. The elements of the latter group g(t) 
are constructed from the elements g of conventional spin rotations in the laboratory frame: 

g{t) = 0-\z, tULt) g 0(Z, LOLt) . (9.5) 

Here the matrix Oai3{z,ut) describes the transformation from the laboratory frame into 
the rotating frame - this is the rotation about the magnetic field axis z by angle uit. 
Now we can find all the degenerate coherent states of the Larmor precession applying the 
symmetry group G to the simplest equilibrium state of the given superfiuid phase: A = 
0~^R*^^''0A*^'^)(R(^))^^, where R^^^ is the arbitrary matrix describing spin rotations in the 
precessing frame and R^^-* is another arbitrary matrix which describes the orbital rotations 
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in the laboratory frame. In case of ^He-B this state corresponds to the state of Cooper pairs 
with L = S = 1 and the total angular momentum J = |24| : 

4°)=Ab5„, . (9.6) 

The action of elements of the group G on this stationary state leads to the following general 
precession of ^He-B with the Larmor frequency (if the spin-orbit interaction is neglected): 

A^i{t) = AR^i{t) , (9.7) 

Rai{t) = 0^p{z,-cut)R^f;}p^^iz,ut){R^%l . (9.8) 
The matrix R*^^-* determines the direction of spin density in the precessing frame: 

5„ = X^i'X ' (9-9) 

where x is the spin susceptibility of ^He-B. This corresponds to the precession of spin with the 
tipping angle cos/9i = R^J'J . The matrix R*^^) determines the direction of orbital momentum 
density in the laboratory frame: 

L, = -R^,{t)S^it) = -xR?jH^ , (9.10) 

(2) 

with the tipping angle cos/32 = Rzz ■ 

C. Spin-orbit interaction as perturbation 

The spin-orbit interaction couples the spin and orbital components of the matrix Ad. 



For ^He-B in (1917]) one obtains [24 



Fd = ^X^l (^Ru{t) - =^X^l(^cose{t) + ^^ , (9.11) 

where is the so called Leggett frequency - the frequency of the longitudinal NMR; 6 is 
the angle of rotation in the parametrization of the matrix Rd in terms of the angle and 
axis of rotation [24]; we shell use here the system of units in which the gyromagnetic ratio 
7 for the '^He atom is 1, hence the magnetic field and the frequency will have same physical 
dimension. 

In the general state of the Larmor precession (19.71) . the spin-orbit interaction contains 
the time independent part and rapidly oscillating terms with frequencies ul, ^col, Sul and 

4 

FD{j)=Fo + Y,FnCos{nuLt) . (9.12) 

n=l 

The time-independent part - the average over fast oscillations - gives the spin-orbit potential 
in fl3.16p . which determines the phases of magnon BEC in '^He-B: 

Fo = Fso(s,/,7) = 

^X^llisl - i + 1(1 + s)(l + /) COS7)' + 

1(1 - sfil - If + (1 - s'){l - + COS7)] . (9.13) 
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Here s = cos /5 and / = cos/3l are z projections of unit vectors s = S/S and 1 = — L/L; 
and 7 is another free parameter of the general precession. Ahogether the free precession is 
characterized by 5 independent parameters coming from two matrices R*-^-* and R*-^^ j4ol |: 
two angles of spin S, two angles of the orbital momentum 1, and the relative rotation of 
matrices by angle 7. In the case of the stationary (non-precessing) magnetization, the 7- 
mode corresponds to the longitudinal NMR mode. 



D. Parametric instability of HPD 

In the simplest description, the dynamics of the 7-mode is determined by the following 
Lagrangian: 



(9.14) 



Here we used the approximation of an isotropic speed of spin waves c. In the time-dependent 
part of Fd we only consider the first harmonic, i.e. according to Eq. fl9.ip we discuss the 
parametric excitation of two 7-modes with ck ^ ul/^. The amplitude of the first harmonic 
is: 



2sl - 1 



Fi = —x^l sin P sin Pl cos(7/2) x 
15 

(1-S)(1-/) , ,^ \ 

^ + (1 + s)(l + /) COS7 . 



(9.15) 



Further we assume that the system is in the minimum of the dipole energy Fq as a function 
of 7. The equilibrium value 7 = 70 is 



COS70 



(2s/ - 1) + 2(1 - s)(l - I) 



(9.16) 



which is valid if the right hand side of Eq. (I9.16P does not exceed unity, i.e. when s + l — 5sl < 
2. 

For the discussion of Suhl instability we need the time-dependent term which is quadratic 
in 7 — 70- Then the Lagrangian 09.141) which describes the parametric instability towards 
decay of Larmor precession to two 7-modes with kc ^ ijJl/2 is (after the shift 7 — 7o — ^ 7; 
neglecting Qi compared to ul', and neglecting the anisotropy of the spin- wave velocity): 



^=2^ (-7' + c'( V7)' + aQh^ cos ULt) 
where, if s + / — 5s/ < 2, the parameter a is 



(9.17) 



a = — sin /3 sin Pl 
15 



3(s + / -s/) 
2(l + s)(l + /) 
35 



1/2 



;i + s)(l + /) + 2(2s/ - 1) + —(1 - s){l - /) 

o 



(9.18) 



39 



Let us rewrite the Lagrangian (19.170 in terms of Hamiltonian as function of creation and 
annihilation operators 6k and 

= /T^-Tm (^'- - ' ^'(^) = ^'k' (9.19) 
= X7k = Vxw.(fc)/2(6k + 6^) , (9.20) 

k 

+ X] (e-^"^*6k&-k + e^'^-*6*61k) , (9-21) 

where we neglected Ql compared to w^,. The spectrum of the excited mode is 
6k (t) = 6ke~^'^^/2*+^''''* , 

//k = ^jMk)-UL/2y-a^Ql/co^,{k) (9.22) 
At the resonance, i.e. when Us{k) = ujl/2, the mode grows exponentially: 

6k(t) oc e^* , A = 2anl/uL ■ (9.23) 

At finite temperatures this growing is damped by dissipation, but at low temperature the 
dissipation becomes small and catastrophic relaxation occurs. Following Ref. |90| one may 
assume the spin diffusion mechanism of dissipation. In this case the equation for temperature 
Teat below which the instability of the homogeneous precession towards radiation of spin 



waves with Us{k) = ck = ujl/2 starts to develop is [9l|, |92[ |: 

D(Teat) = 2\cycul . (9.24) 

Here D{T) is the spin diffusion coefficient, which depends on temperature and decreases 
with decreasing T. This agrees with observations. 

X. BEYOND MAGNON BEG: SPIN SUPERGURRENTS AND SPIN HALL EF- 
FEGTS 

A. Microscopic theory of spin supercurrent in ^He-B 

Here we give the 'microscopic' derivation for the spin supercurrent, which has been dis- 
cussed on the phenomenological level of magnon BEC. The underlying microscopic physics 
is the BCS theory of p-wave spin-triplet superfiuid '^He. Superfiuid spin currents in ^He-B 
exist even in the absence of magnon BEC. They come from the spontaneous breaking of 
spin-rotation symmetry. The spin supercurrent in ^He-B is expressed in terms of the spin 
superfiuid velocities: 

^ai = '^ida = -^eafSjR/Bj^ iRjj ■ (10-1) 

The corresponding gradient energy is 

Tgrad = -^Pai,lij^ai^Pj , (10-2) 
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where Pai,i3j is the spin rigidity tensor with spin rigidity parameters 

Pai,l3j = ^ [Cy^Q^^jj — (c| — c\){RaiRi3j + RajR(3i)] ■ (10.3) 

From these equations one obtains the spin supercurrent in ^He-B: 

Jai = 77^ = —Pai,l3j^l3j ■ (10-4) 

This spin current averaged over the fast precession determines the parameters of the phe- 
nomenological equations f l4.16p and fl4.17p for the spin supecurrent emerging in magnon 
BEC. 



B. Spin current and electric field in '^He-B 



The symmetry properties of the spin superfluid velocity, and thus of the spin supercurrent, 
allows to couple linearly the spin current with electric field even in t he absence of the spin- 
orbital interaction. The following term in the action is possible |l01 |: 



F /3cijj^ujQ,iRQ,j El 



(10.5) 



The parameter /3 is not well defined from the microscopic theory due to t he unknown Fermi- 
liquid corrections involved. The estimate for /3 reported in Ref. 102| is /3 ~ 10~^e/cm. 
Variation with respect to Uai demonstrates that there exists a linear response of the spin 
supercurrent on the electric field: 



J rv 



dF 



/Scijf^Rf^jEjf 



(10.6) 



This spin current is transverse to the electric field and thus represents the spin current Hall 
effect. As distinct from the spin Hall effect predicted by Dyakonov and Perel 
spin-Hall effect does occurs in the absence of spin-orbit interaction. 



103 , this 



C. Electric and magnetic fields as SU{2) gauge fields 



The interaction of the electric and magnetic fields with the order parameter in superfluid 
phases of '^He, can be also found using observation that H and E may be considered as 
temporal and spatial components of the SU{2) gauge field, where SU{2) is the group of 
the spin rotations. The auxiliary SU{2) gauge field A'^ is convenient for the description of 
the effects related to the spin current. In the spinor representati on one ha s the following 



covariant derivatives coming from the auxiliary SU{2) gauge field [lOll . |104 



D = V - iA* 



D^ = dt + i^^ 



cr- 



(10.7) 



Some components of the field A'^ are physical, being represented by the real physical quanti- 
ties which couple to the fermionic charges. Example is provided by the Pauli magnetic field 
H\ which play the role of the component Aq of the SU{2) gauge field, see ( 19. 3 p jiol . 105 
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while the spatial components are pl ayed by the ele ctric field which enters the gradient energy 
via the covariant spatial derivative (lOl . 104 . 106 1: 



A!. 



Eh . 



(10.8) 



The electric field enters due to the relativistic spin-orbit interaction of the spin of ^He atom 
with the electric field E. 

The spin current is obtained as variation of the action with respect to the fictituous SU (2) 
gauge field: 



5S 



10.9) 



After the spin current is calculated the values of the auxiliary fields are made equal to zero 
or to the values of the corresponding physical fields which simulate the gauge fields. 
For example, in the presence of electric field, equation fll0.2l) becomes 



1 1 

—e^aikEkji^Pj — —ei3jiEi) 



(10.10) 



which demonstrates that electric field enters as the SU (2) gauge field forming the covariant 
derivative. This gives another response of the spin cupercurrent on the external electric 
field: 



dK 



grad 



7 



dcOn 



~Pai,l3j{^/5j ^IBjkEk) , 



;io.ii) 



As distinct f rom (110. 6p this spin-Hall effect is governed by the spin-orbit interaction. Accord- 
ing to Ref. 1071] both spin Hall effects in ^He-B should modify the spin current Josephson 
effect in magnon BEC. The supercurrent, induced by electric field, leads to an additional 
phase shift proportional to electric field, which is to be measured. 



D. Quantum spin Hall effect 

There are several types of responses of spin and electric currents to transverse forces 
which are quantized in 2+1 systems under appropriate conditions. The most familiar is the 
conventional quantum Hall effect (QHE). It is quantized response of the particle current to 
the transverse force, say to transverse gradient of chemical potential, J = a^yZ x V/i. In 
the electrically charged systems this is the quantized response of the electric current to 
transverse electric field J*^ = e'^a^yZ x E. 

The other effects involve the spin degrees of freedom. An example is the mixed spin 
quantum Hall effect: quantized response of the particle current J (or electric current 3^) 
to transv erse gradient of magnetic field interacting with Pauli spins (Pauli field in short) 



108, 109 



J = cr™z X V{-fH') , J-^ = eJ . (10.12) 

The related effect, which is determined by the same quantized parameter crj!^™'^, is the 
quantized response of the spin current, say the current of the z component of spin, to the 
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gradient of chemical potential |110l ]. In the electrically charged systems this corresponds to 
the quantized response of the spin current to transverse electric field: 



^mixcd^ X V/i 



xy 



mixed 
xy 



Z X E, 



;io.i3) 



This kind of mixed Hall effect is now used in spintronics . 

Ffinally there is a pure spin Hall effec t - the quantiz ed response of the spin current to 



transverse gradient of magnetic field |l08l . Il09l . Ill2l . |113 



spin/spin 
xy 



z X V(7/7" 



(10.14) 



Let us consider the mixed spin Hall effects in f llU.13p and f llU.14p . These two effects are 



related, since they described by the same topological Chern-Simons action |108j and thus 



by the same parameter o"^^^^^. To see this, let us remind that the spin current is obtained 

For 



as variation of the action over the fictituous SU{2) or 5*0(3) gauge field, see (llO.Qp 
example, the current of the z-projection of spin is 



6S 



The corresponding Chern-Simons term in the action is given by [108 



;io.i5) 



d^xdtAlVryA 



(10.16) 



where Ap is the vector potential of the conventional electromagnetic field, and Al = {Af, Aq) 
represent components of auxiliary (fictituous) SU{2) gauge field. Variation of the action with 
respect to the field Al gives the spin current in (110. 13p . On the other hand, the variation 
of the action with respect to field Ai gives electric current in terms of the gradients of an 
auxiliary gauge field. However, from equation (19. 3p or (110. 8p it follows that the role of the 
auxiliary gauge field A^ is played by magnetic field ff^. As a result one obtains equation 
(110. 12p for electric current. 

Equation (110.160 ha s be en originally introduced for a thin film of the so-called planar 



phase of superfluid He |108| . However, it is better suited for the two-dimension al to pological 
insulators with time reversal invariance (on topological insulators see review |ll4l |). These 
materials have the same topological structure as the planar phase, which is also time reversal 
invariant, but the advantage of these materials is that the are insulating and thus the 
superconductivity does not mask the spin Hall quantization. 

Discussion of the mixed Chern-Simons term can be found in Ref. 



115 . For the related 



phenomenon of axial anomaly in particle physics, the mi xed a ction in terms of different (real 



and fictituous) gauge fields has been introduced in Ref. [116 



XI. CONCLUSION 



The superfluid phases of liquid ^He at extreme low temperatures are unique states of con- 
densed matter with physical properties which can be compared to the vacuum of relativistic 
quantum field theories. The A phase (^He-A) belongs to the same symmetry and topology 
class as the vacuum of the Standard Model of particle physics in its massless {i.e gapless) 
phase and can also be described as a semi-metal-like system with non-trivial topology. The 
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B phase (^He-B), in contrast, is similar to the vacuum of the Standard Model in its massive 
(or gapped) phase and to 3-dimensional topological insulators with time reversal symmetry. 
In addition to fermionic quasiparticle excitations, superfluid ^He has also bosonic quasipar- 
ticles, such as magnons - quanta of excitations of the magnetic subsystem. These magnon 
excitations can form long-lived Bose-Einstein condensates both in ^He-A and ^He-B, and 
these condensates experience their own superfluidity, which is not related to superfluidity of 
the underlying system. 

Formally, the phenomenon of superfluidity requires the conservation of charge or particle 
number. However, the consideration can be extended to systems with a weakly violated 
conservation law, including a system of sufficiently long-hved quasiparticles - discrete quanta 
of energy that can be treated as real particles in condensed matter. The spin superfluidity 
- superfluidity in the magnetic subsystem of a condensed matter - is manifested as the 
spontaneous phase-coherent precession of spins first discovered in 1984 js], llf. This superfiuid 
current of spins is one more representative of superfiuid currents known or discussed in other 
systems, such as the superfiuid current of mass and atoms in superfiuid ^He; superfiuid 
current of electric charge in superconductors; superfiuid current of hypercharge in Standard 
Model; superfiuid baryonic current and current of chiral charge in quark matter; etc. The 
analogy of the dynamical superfiuid state of coherent precession with the non-perturbative 
dynamics of the physical vacuum is to be found. 

Different condensates and thus different states of magnon superfluidity have been created 
by choosing different experimental arrangements. At low temperatures the condensate is 
conflned in a magnetic trap which is formed by the order parameter texture of the superfluid 
state. This produces the analog of atomic BEC in laser traps, but adds some new features, 
such as formation of the non-ground-state condensate and the mass of the Goldstone bosons. 
The magnon condensates can be used to probe the quantum vacuum of '^He in the limit 
r — )■ 0, where conventional measuring signals become insensitive. 
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